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Abstract 

We compute asymptotics for Hankel determinants and orthogonal polynomials with respect to a discon¬ 
tinuous Gaussian weight, in a critical regime where the discontinuity is close to the edge of the associated 
equilibrium measure support. Their behavior is described in terms of the Ablowitz-Segur family of solutions 
to the Painleve II equation. Our results complement the ones in [33]. As consequences of our results, we 
conjecture asymptotics for an Airy kernel Fredholm determinant and total integral identities for Painleve II 
transcendents, and we also prove a new result on the poles of the Ablowitz-Segur solutions to the Painleve II 
equation. We also highlight applications of our results in random matrix theory. 


1 Introduction 


Consider the Hankel determinant, 


= detj / x^*^w{x)dix 
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with respeet to a diseontinuous Gaussian weight of the form 
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The weight is periodie in jS and we ean restriet to the ease -1/2 < Re^S < 1/2 without loss of generality. If [A is 
purely imaginary, the weight is positive. 


We also eonsider the monie orthogonal polynomials p„ of degree n with respeet to the weight w(x) on the real 
line, defined by the orthogonality conditions 


Pnix)Pmix)w(x)dx = h 

n^n 


hn P). 


(1.3) 
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1 Introduction 


Those are conneeted to the Hankel determinant by the well-known identity H„{Aq,P) = OLo hk(Ao,/3). We 
denote by R„ = R„(Ao,/3) and Q„ = Qn(Ao,/3) the recurrence coefficients in the three-term recurrence relation 

XPnix) = p„+l(x) + Q„Pn{x) + RnPn-\{x). (1.4) 

The question which we are concerned with in this paper is the large n behavior of the Hankel determinants 
Hn, the polynomials Pn(x), and their recurrence coefficients and Q„, in the regime where the point of 
discontinuity To behaves like V^. They can asymptotically be expressed in terms of the Ablowitz-Segur 
solutions to the Painleve II equation. As important by-products of the asymptoties for the Hankel determinants, 
we also conjeeture so-ealled large gap asymptotics for an Airy kernel Fredholm determinant and total integral 
identities for the Ablowitz-Segur solutions of the Painleve II equation. Relying on a result of [33], we in 
addition prove a new result about the poles for those Painleve transcendents. 

If we let To = T V^, the large n asymptotics of the orthogonal polynomials (1.3), the recurrence coefficients 
(1.4), and the Hankel determinant (1.1) depend dramatieally on whether |T| < 1 or |T| > 1, i.e. whether the jump 
loeation To is inside or outside of the support [- ^/2n, V^] of the equilibrium measure with Gaussian external 
field. In the case |T| > 1, all the objects of interest behave effectively as they do for the pure Gaussian weight 
(i.e., the case where we formally set To = -i-oo); the discontinuity yields an exponentially small eorreetion only 
[25]. In the case |T| < 1, the situation is different; the discontinuity of the weight becomes strongly visible in 
the large n behavior of the orthogonal polynomials, the recurrence coefficients, and the Hankel determinant 
[21]. For the Hankel determinant, it was proved in [21, equation (1.5)] that 

HniAo,/3) = Hn(Ao,0) G(1 +/3)Ga-m - x 

X exp(2myS(arcsinT-r TVr^^))|lIReySI < ^, (1.5) 

as n ^ oo, uniformly for T in compact subsets of (-1,1). Here G is the Barnes’ G-function, and 


n-\ 

HAAo,0) = i2nf^2-""/^Y]^- 

k=i 


( 1 . 6 ) 


2 

denotes the Hankel determinant eorresponding to the pure Gaussian weight e~^ . Asymptoties for the reeurrence 
coefficients Qn and in the case -1 < T < 1 are also given in [21]. 

In this paper, we analyze the transition regime where the point To of discontinuity of the weight is (relatively) 
close to V^. More precisely we let 


To = T V^, 


T = 1 -I —n 
2 


-2/3 


(1.7) 


where r 6 R. We will see that the asymptotic behavior of H„, pn,Rn, and Qn depends in a non-trivial way on 
the parameter t in (1.7). The asymptotic behavior is described in terms of a family of solutions to the Painleve 
II equation 

Utt = tu + 2u^, ( 1 . 8 ) 


with the asymptotic behavior 


u{t\ k) ~ /fAi(r), t -1-00, 


(1.9) 


where Ai denotes the Airy function, and 


u(t; k) 


1 


yj2ifi sin ^(t; jS) +0 


1 


^2-3|Re/?| I’ 


( 1 . 10 ) 
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with 


- ilog - ^z^log(-t)-3z^log2, zc^ = 1 - e |RejS|<^. (1.11) 

For 0 < zc < 1, these solutions are known as the Ablowitz-Segur solutions [1] of the second Painleve 
equation. They are uniquely eharaeterized either by (1.9) or by (1.10). Moreover, it is known that zz(t; zc) has 
no singularities for t on the real line if zc 6 zR or if |zc| < 1. For zc 6 R \ [-1,1], or equivalently | RejS| = 1 /2, 
it is known that zz(r; zc) does have real poles [4]. Relying on a result from [33], we will prove the following 
result, stating that u has no real poles for any zc 6 C \ ((-oo, -1] u [1, +oo)), or equivalently for any jS with 
\Rej3\ < 1/2. 

Theorem 1. Let u(t;K) be the solution to the Painleve II equation (1.8) characterized by (1.9). If k e C\ 
((-oo, -1] u [1, +oo)) is fixed, then u{t; zc) has no poles at real values oft. 

In the ease zc = 0, we simply have zz(t; zc) = 0; the unique Painleve II solution satisfying (1.9) with zc = +1 
(whieh means formally that jS = -zoo) is known as the Hastings-MeLeod solution. 

The funetion y(t; /3) = zz(t; zc)^ solves the Painleve XXXIV equation 

7 (y*)^ 

ytt = ^y +'lty+f—. (1.12) 

2y 

The funetion y{t,/3) and equation (1.12) are, in faet, the objeets whieh direetly appear in our double sealing 
analysis of //„, and Q„. Our next result deseribes the asymptoties or the Hankel determinants H„{Aq,/3). 

Theorem 2. Let | Re^SI < 1/2 and let H„{Ao,j3) be the Hankel determinant (LI) corresponding to the weight 
(1.2), with Ao given by (1.7). If = 1 - we have 


//„(do,;0) = e'"^”//„(do,O)exp 


(T-t)u(T\K) dr (1 +o(l)). 


oo. 


(1.13) 


\ j t / 

uniformly for t 6 [-M, co) for any M > 0 and for fi in compact subsets o/| Re;S| < 1/2, where //„ (Aq, 0) is 
given in ( 1 . 6 ). 


Theorem 2 has two eonsequenees whieh are not directly related to the Hankel determinants or orthogonal 
polynomials studied in this paper, but whieh are of independent interest. To deseribe them, we note first 
that the exponential in (1.13) can be reeognized as the Traey-Widom formula for the Fredholm determinant 
det^l - where is the integral operator with kernel 


^Ai(vy) 


Ai(x)Ai'(y)-Ai(y)Ai'(x) 
x-y 


(1.14) 


aeting on \t, +oo). Indeed, it was shown in [31] that 


det(l 


■ 




(t - t)u(T-, KydiT 


(1.15) 


This observation, together with a strengthened version of the Hankel determinant asymptoties (1.5), allows us 
to formulate the following eonjecture about the t -oo asymptotics of det^l - zc^XaiIj^^^^). 
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Conjecture 3. Let /c 6 C \ ((-oo, -1] u [1, +oo)) and define /3 by 


1 - 


\Re/3\ < 1/4. 


(1.16) 


As t ^ -oo, we have 

logdet(l - - ^>e"log(-0 + log(G(l +fi)G{\ -(5)) - 3ldHogl + o(l), (1.17) 

or equivalently inform of a total integral identity 

Hm {T-t)u{T-,KfdT + + ^yS^log(-t)j = log(G(l +y8)G(l -yS)) - 3fHog2. (1.18) 

Remark. Similar asymptotics for the Airy kernel determinant in the case k = \ were proved in [2, 12]: we then 
have 

fi 1 

log det(/-=-pr - ^ logit] + Co+ G(r^), t^-oo, (1.19) 


i[f,+oo)/ 12 8 

where co = log 2/24 + ^'(-1) and ^ is the Riemann ^-function. As /r ^ 1, it was shown recently in [7] that 

12 8 

as long as /c < 1 , and k ^ I sufficiently rapidly so that 


f 1 

log det (/ - - log |t| + Co + o ( 1 ) ^ 


t -OO, 


( 1 . 20 ) 


log(l - K^) 2 V 2 

(-t)3/2 ^ ^ 


( 1 . 21 ) 


The total integrals of different expressions involving the second Painleve transcendent were studied in [3] . The 
integral (1.18) does not belong to the type which can be handled by the technique of [3] . Indeed, like the similar 
integral corresponding to equation (1.19), the integral in (1.18) belongs to the third, most difficult type of total 
integrals of Painleve functions as classified in the end of Section 6 of [3]. This means that the evaluation of 
this integral goes beyond the analysis of the Riemann-Hilbert problem corresponding to the Ablowitz-Segur 
Painleve II transcendent. As we already indicated, the proof of (1.18) can be achieved via an improvement of 
the error term in (1.5). Another possibility is to use certain differential identities for the Airy determinant in 
(1.15) with respect to k. We intend to consider these issues in our next publication. 

Additionally, the asymptotics of the PXXXIV transcendent y (t;yS) = u(t;K) as t ^ -oo can be calculated 
directly by the same method as the ones for t +oo. We will not present this computation since it is mostly 
identical to the one in [21], and alternatively this asymptotics can be obtained using the connection formulae 
for the Painleve II equation [26]. Moreover, the following singular asymptotics take place when Re/3 = 1/2. 

Theorem 4. Let u(t; k) be the solution to the Painleve II equation (1.8) characterized by (1.9) and let = 
1 - = 1 + j3 = \/2 + iy, y e R. Then y (t; fi) = u (t; k)^ is a solution to the Painleve XXXIV equation 

(1.12) and has the following asymptotics as t ^ -oo, away from the zeros of trigonometric functions appearing 
in the denominators: 


y\L:: + iy\ = 


-t 1 

+ 


cos2 0 ^f22t 


1 _ 2 y 3(1272 - i)sin0 

-7 + - tg 0 + -^- — - 

2 cos2 0 16cos^0 


+ 0 


where 0 (t; 7 ) = - {-tf^^ + ^ 7 log (-0 + 37 log 2 - arg TI ^ + f 7 


1 


( 1 . 22 ) 
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Asymptotics of this type in relation to the seeond Painleve equation have been previously obtained via different 
methods in [26] and in [9], but the second term is a new result of the present work. This computation is based 
on an undressing procedure adapted from [9]. Thus we will not present the derivation of (1.22) either. 

For the recurrence coefficients and Q„, we have the following result, which was partially obtained before in 
[33], see Remark 7 below. 

Theorem 5. Let Rn and Qn be the recurrence coefficients defined in (1.4), associated to the orthogonal 
polynomials with respect to the weight ( 1.2). Let | Re/3| < 1/2 and let Aq be given by (1.7). Then, as n ^ oo, 
the recurrence coefficients have the following expansions, 

Rn(Ao,/3) = ^ - ^u(f,K)^n^/^+0(l), (1.24) 

and 

Qn(A^ffi) = -\u{f,Kfn-^l^+0(n-^l^), (1.25) 

V2 ^ ’ 

uniformly for t 6 [-M, oo]for any M > 0, where k is given by (\.16). Additionally, we have the asymptotics of 
the normalizing coefficients hr,: 

+n~^^^v(t;K)dt + n~^^^^(v(t;K)'^ - m (t;/^)^) + , (1.26) 

where 

/ oo 

u (r; k)^ dr. 

Remark 6. The formal substitution, 

t =-2(1 - A)n^'^ 

in the asymptotics for the recurrence coefficients transforms them, with the help of the asymptotic expansion 
(1.10), into the non-critical asymptotics obtained in [21]. This important fact indicates, at least on the formal 
level, that the description of the transition regime in the large n behavior of the recurrence coefficients is 
complete. 

Remark 1. The general form of (1.24) and (1.25) was formally suggested in [23] (together with the asymptotic 
characterization of the Painleve II function u(t-, k)) and it was proved by Xu and Zhao in [33]. They obtained 
their asymptotic expansions in terms of u(t) = 2^^^m(2“'^^t)^. It was noted that this is a solution of a Painleve 
XXXIV equation, but no asymptotics for u(t) as t ^ +00 were obtained, and thus the authors of [33] did 
not identify u in terms of the Ablowitz-Segur solution characterized by (1.9) or (1.10). In fact, assuming the 
matching of the estimates (1.24) and (1.25) with the non-critical formulae of [21], asymptotics for u(t) as 
t -00 were deduced heuristically. There is, however, no independent derivation of it which is needed for 
the rigorous completion of the analysis of the transition regime in question. The -l-oo - characterization of the 
Painleve transcendent u(t), even heuristically, is not given in [33]. 

As an additional result, we also obtain an analog of the Plancherel-Rotach asymptotics for classical Hermite 
polynomials [29]. 

Theorem 8. Let Pn(x) be the degree n monic orthogonal polynomial with respect to the weight (1.2), and let Aq 
be given by (1.7). Let |Re/3| < 1 /2. Then, as n ^ 00 , 

Pn(do) = -^(y) n^'^e''^'^\(t\K)(l+0(n~^''^)), (1.27) 

with K given by (\.\6). 
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Remark 9. Using the asymptotic behavior (1.9) for u as t ^ +oo, (1.27) matches formally with the classical 
Plancherel-Rotach asymptotics for the Hermite polynomials [29]: 

Pn(Ao) = (0 (i n ^ oo, (1.28) 

where p„ are the monic Hermite polynomials. 

On the other hand, if we let yS ^ 0, or equivalently /r ^ 0, we have (see equations (5.29) and (5.31) below) that 

u(t;K) = 0, lim-M(t;A:) = Ai(t), (1.29) 

K 

and this allows us to recover (1.28) also in this limit. 

Remark 10. Consider the case of purely imaginary yS, i.e. yS = iy, y 6 R. Then the three-term relation (1.4) 
generates in the usual way (see e.g. [11]) a symmetric on I 2 Jacobi operator, LP, defined by the semi-infinite 
matrix 

= R]!+\Sn+\,m + Qn6n,m + H, m > 0, Rq = 0 

whose domain is D = = {pQ,p\,.. y e I 2 : Pk = 0 for sufficiently large fcj. Since the moment problem for 

the measure 

dyu (x) = w (x) dx 

with the density w{x) given by (1.2) is determinate, the operator LP is essentially self-adjoint and dp(x) is the 
spectral measure of its closure L = LP. Therefore, the results of our last two theorems provide an insight into 
the properties of semi-infinite Jacobi matrices, i.e. discrete Schrodinger operators on a half-line, whose spectral 
densities have discontinuities. In the earlier works [10, 21, 25] it was demonstrated that the discontinuities 
in the spectral density are responsible for the oscillatory pattern in the large n asymptotics of the entries 
of the Jacobi matrix (in the coordinate asymptotics of the potentials of the discrete Schrodinger operator). 
More precisely, the oscillations occur when the point of the jump of the density is inside the support of the 
corresponding equilibrium measure. If the jump is outside, the behavior of the potentials R„ and is monotone. 
Formulae (1.24), (1.25) and (1.27) describe the corresponding transition regime. The formulae show that if 
the jump happens near the edge of the support then the large n (coordinate) asymptotics is governed by the 
Ablowitz-Segur solution of the second Painleve equation and the parameters of the solution are explicitly 
related to the size of the jump. We actually believe that this fact is universal, i.e. the transitional formulae will 
be the same even if the Gaussian background in the spectral measure is replaced by an arbitrary exponential 
weight. 

Our proofs of Theorem 5 and Theorem 8 are based on the nonlinear steepest descent method of Deift and Zhou 
(or, rather on its adaptation [15] to the Riemann-Hilbert (RH) problems related to the orthogonal polynomials 
[19]). This method was applied in [33] to the case of a discontinuous Gaussian weight with the point of 
discontinuity scaled as in (1.7). We will rely on the transformations and results from this paper, but we will 
adapt them in such a way that we can identify the function u{t; k) as the Painleve II solution with asymptotics 
(1.9) and (1.10). The RH analysis is presented in Section 3, and the proofs of Theorem 5 and Theorem 8 are 
given in Section 5. 

Theorem 2 can be proved in two different ways. The first one is very short and relies on the Tracy-Widom 
formula (1.15) and on known asymptotic results in the Gaussian Unitary Ensemble. This proof will be given in 
Section 2. The second proof, given in Section 6, is lengthy but has the advantage of being self-contained. It 
relies on the RH analysis which we need anyways for the asymptotics of the orthogonal polynomials and their 
recurrence coefficients. As is always the case in the asymptotic analysis of Hankel and Toeplitz determinants. 
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the move from the asymptoties for the orthogonal polynomials and its recurrenee eoefficients to the asymptoties 
for the Hankel determinants is nontrivial. One has to address the “constant of integration problem” (c.f. [17]) 
which we do with the help of relevant differential identities for the Hankel determinant Hn{Ao,/3). 

In the RH analysis, we will identify the function u(t; k) as the solution to the Painleve II equation with 
asymptotics (I.9)-(I.I0) using Lax pair arguments and an asymptotic analysis for a certain model RH problem 
(see Section 4), which is equivalent to the one which appeared in [33]. Solvability of this model RH problem 
was proved in [33], and we prove Theorem 1 as a consequence of this in Section 3.5. 

The analysis in this paper shows similarities with the work [22] where a Painleve XXXIV function appeared in 
a parametrix for a different type of critical edge behavior in unitary random matrix ensembles, namely with 
a root singularity instead of the jump singularity which we consider here. The RH problem which we study 
differs, however, from the one analyzed in [22]. This yields, in particular, serious technical differences in the 
analysis of the large positive t behavior of the Painleve transcendent. 

Remark 11. As it has already been indicated, it is Painleve XXXIV equation (1.12) and the corresponding 
model RH problem that appear naturally during the asymptotic analysis of the orthogonal polynomials Pn(x). 
The solution y(t; P) which emerges in this analysis is characterized by its RH data. We need to transform this 
characterization into the asymptotic behavior of y{t',P) as t —> +oo. Because of the relation y = between the 
solutions of Painleve XXXIV equation (1.12) and the solutions of Painleve II equation (1.8), one could think 
that the needed asymptotics could be extracted from the work of A. Kapaev [26], where the complete list of the 
global asymptotics of the second Painleve transcendent is presented. However, to be able to use the results of 
[26] one needs to connect the RH data of y(t) with the RH data of u(t). A well-known though still striking fact 
(see e.g. Chapter 5 of [20]) is that there is no simple relation between the Lax pair and the RH problem for the 
Painleve XXXIV equation (1.12) and the standard Flaschka-Newell Lax pair (which is used in [26]) and the 
RH problem for the Painleve II equation (1.8). Henee one does not know a priori the asymptoties of u{t). There 
exists, however, a simple relation between the Lax pair and the RH problem for the Painleve XXXIV equation 
(1.12) and the Lax pair and the RH problem for the nonuniform second Painleve equation 

q„ = tq + 2q^-^, (1.30) 

so that one can use [26] and determine the asymptoties of q(t). Unfortunately, now the problem with translation 
of the asymptotics of q(t) into the asymptotics of y(t) arises. The fact of the matter is that the relation between 
the Painleve functions y(t) and q(t) is more complicated than the relation between the Painleve functions y{t) 
and u(t). Indeed, one has that 


yit) = U{t) = q\t) + q\t) + ^ (1.31) 

(see e.g. [22, Appendix A]). This formula virtually destroys the asymptotie information which one eould obtain 
for the function q{t) from [26]. For instance, one finds from [26] that the function q{t) behaves as ~ yl-t/2 
as t ^ -oo. This, as we know a posteriori, must translate to the exponentially decaying asymptotics of 
yit) as t ^ - 1 - 00 . It is extremely difficult to verify this directly using (1.31): one has to prove cancellation of 
an asymptotic series in all orders of magnitude. Even worse is the situation with the asymptotics of qit) as 
t - 1 - 00 . It is singular (and is described in terms of the cotangent function) and, after substitution into (1.31) it 
should transform into an oscillatory smooth decaying asymptotics. We refer the reader to [22], where a similar 
phenomenon had already been encountered, for more details. The above discussion makes it clear that, in 
spite of the simple relation to the second Painleve function m( 0, an independent asymptotic analysis of the 
Painleve XXXIV funetion yit) is necessary. Of course, it is enough to evaluate the asymptotics of y(L/?) either 
for t - 1-00 or for t -oo, since the one-end asymptotics will enable us to identify the function uix) and use 
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[26] to determine its asymptoties on the another end. We have chosen to evaluate the asymptotics of y(t;/3) as 
t + 00 . The relevant nonlinear steepest descent analysis is presented in Section 4. This analysis has some 
new teehnieal features whieh are specifically indicated at the beginning of Section 4. 

1.1 Applications 

We eonelude this introduction by indicating some applieations of our results. 


1.1.1 Random matrix moment generating function 


Consider the n-dimensional Gaussian Unitary Ensemble (GUE) normalized such that the joint eigenvalue 
probability distribution is given by 

1 " 2 

— Y\ ~ FI a:i,..., G R. (1-32) 

” \<,i<j<Sn 7=1 

The partition funetion Z„ is then equal to «!//„ (do, 0), with H„ (do, 0) given in (1.6). For an nxn GUE matrix, 
define the random variable as 

= number of eigenvalues greater than do. (1.33) 


It is natural to ask how the average of „ or its variance behaves for large n. The Hankel determinant with 
a discontinuous Gaussian weight carries information about such quantities. Indeed, the moment generating 
function of the random variable whieh is defined as M^g_„(y) := E„ " j, can be expressed as 


MAQ,niy) 


1 

z„ 



dn 

7=1 


2 X 



Xj < do 
Xj > do 


'V 

X dXy 

/ 


(1.34) 


This is in faet the ratio of two Hankel determinants, one with a discontinuous Gaussian weight, and one with a 
regular Gaussian weight: if we write y = -2ni/3, we have 


MAo,n(y) 


Hn(Ao,0) 


(1.35) 


This is true for any n and do. 

The large n asymptoties for the Hankel determinant proved in Theorem 2 together with the explieit expression 
(1.6) for Hn (do, 0 ), immediately give information about the moment generating function as n ^ oo if do is 
scaled as in (1.7). 

Expanding the moment generating funetion for small values of y, we have 

MA,,niy) = 1 +yE„(X,„,„) + ^E„(X2^ J +G(y3), y ^ 0, (1.36) 

so the average and variance of „ can be read off immediately from the small asymptotics for the Hankel 
determinant. 

In particular, differentiating (1.13) with respeet to /3 and using (1.35) and (1.36), we obtain 



1 d^ / 


(r - t)u{T\ 


/J=0 


(1.37) 
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with K given by (1.16), which means that the large n limit of the moments of the random variable „ can 
be expressed in terms of the Ablowitz-Segur Painleve II solutions u(t; k) and its /r-derivatives evaluated at 
/<■ = 0. Note that this differentiation is justihed since the asymptotics (1.13) are known to be uniform in a small 
neighborhood of /3 = 0. The hrst /c-derivative of u is the Airy function, by (1.29), and this implies that 


lim E„(Xio,„) 


/ +00 ^ 

(r - OAi (T)2dT = - {tf - Ai (r)Ai '{t) - 2tA\ '{tf ). 


(1.38) 


The same formula can be derived directly from lim„.^ooE„ (X^o,„) = p{T)dT , where pit) = ^Ai(b0 = 
Ai '{t)^ - Ai "(OAi it) is the density for the largest eigenvalue. Similarly, the behavior of higher moments can 
also be studied via just the correlation functions p„, (vi,..., x„,) = det (^K^i (v,, Vyjj, We would like to thank 
Peter Forrester for pointing out this fact. 


1.1.2 Largest eigenvalue in a thinned GUE 

The second application is connected to the so-called thinning procedure in the GUE. Consider the n eigenvalues 
> ... > Xn of a GUE matrix, and apply the following thinning or filtering procedure to them: for each 
eigenvalue independently, we remove it with probability s 6 (0,1). This leads us to a particle configuration, 
where the number of remaining particles can be any integer i between 0 and n, and we denote those particles 
by pi > ... > pf. Below, we show that the largest particle distribution in this process can be expressed in terms 
of a Hankel determinant with discontinuous Gaussian weight. More precisely, we have 


Prob, (pi < To) = s). 


(1.39) 


where M^^^nit) is defined in (1.35). 

To prove (1.39), write E„ik, To) for the probability that anxn GUE matrix has exactly k eigenvalues bigger 
than To. If we want none of the thinned or filtered particles pi,... ,p^ to be bigger than To, that means that all 
GUE eigenvalues which are bigger than To have to be removed by the thinning procedure. Therefore, we have 

n 

Prob,(pi < To) = Enik, To)/, (1.40) 

k=0 

since each eigenvalue is removed independently with probability s. 

Using the integral representation (1.34), it is on the other hand straightforward to show that 

n 

Mi0.n(log/ = Xi^/^’^0)/. (1.41) 

k=0 


Alternatively, this follows from the equation 


Enik, To) 


1 

U 



M^o,„(log s). 


(1.42) 


which is well-known and proved, for example, in [28, Ch. 6 and 24]. Combining (1.40) with (1.41), we obtain 
(1.39). Consequently, by (1.35) and (1.13), 

lim Probi(pi < To) = lim „(log s) = exp 


(r - t)uiT; K)^dT 


(1.43) 
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2 Theorem 2 and Conjeeture 3 


where s = \ - k^. This relation, without the Hankel determinant, was discussed previously in [5, 6], where a 
transition was observed from the Tracy-Widom distribution (at 5 = 0) to the Weibull distribution (at 5=1). 
It is challenging, however, to describe explicitly the transition asymptotic regime from the behavior (1.17) 
corresponding tOyS = -^ In 5, 0 < 5 < 1 to the Tracy-Widom asymptotic behavior, 

lndet(l-i^Ai|[^^„)) = ^?^-^ln(-0 + ^ln2 + ^'(-l) + o(l), t ^-oo, (1.44) 

corresponding to 5 = 0, i.e. p = -ioo ox k = 1. Here, ^ is the Riemann zeta-function. Similar transition 
regime for the sine-kernel determinant has been already described in [8] in terms of elliptic functions, and the 
presence of a very interesting cascade-type asymptotic behavior has been detected (see also [18] where the 
problem was analyzed, on a heuristic level, for the first time ). In the case of the Airy-kernel, the question is 
still open, although on the level of the logarithmic derivatives, i.e. on the level of the Painleve function u {t; k), 
the transition asymptotics from the Ablowitz-Segur case {k < 1) to the Hastings-McLeod (/c = 1) case has 
already been found in [7]. 

1.1.3 Random partitions 

The Airy kernel Fredholm determinant can be interpreted in terms of random partitions. The Plancherel 
measure on the set of partitions of A 6 N is a well-known measure which has its origin in representation theory. 
It can be defined in an elementary way by the following procedure. Take a permutation cr in 5 ^ and define 
Xi as the maximal length of an increasing subsequence of cr. Next, we define X 2 by requiring that xi + X 2 is 
the maximal total length of two disjoint increasing subsequences of cr. We proceed in this way, and define Xk 
recursively by imposing that -l- ■ • • -l- Vjt is the maximal total length of k disjoint increasing subsequences 

of cr, and we continue until Xi -\ - + xt = N. This procedure associates a partition xi >■■■> x„ of N to a 

permutation cr e S n. The uniform measure on 5 ^ induces a measure on the set of partitions of N, which is the 
Plancherel measure. 

We now take a random partition xi > • ■ • > x„ of N with respect to the Plancherel measure. Then, the particles 
N~^^^(xi - 2 ^fN) converge to the Airy process as N ^ oo, see e.g. [30]. Therefore, if we apply the filtering 
procedure which removes each component x, of the partition independently with probability s, we obtain a new 
partition ni > ■ ■ ■ > ix,„ of a number C < A. Using similar arguments as in [30], it follows that 

Prob,(A-i/^0ui - 2 VA) < t) = det(l - (1 - ^) ■ (1-45) 

Note that the conjectured integral identity (1.18) is of value in relation to (1.43) and (1.45). 

We want to conclude this section by making the following general remark. From the point of view of the random 
matrix theory the examples considered in this section indicate that, in fact, it is the whole Ablowitz-Segur 
family of the Painleve II transcendents that could appear in the theory and not only the Hastings-McLeod 
solution. Regarding the second appearance, it has already been known due to Bohigas et. al. [5, 6], however 
the first and the third examples are apparently new. 


2 Theorem 2 and Conjecture 3 

2.1 Proof of Theorem 2 

Denote for the GUE eigenvalue correlation kernel 

n-l 

KAx,y) = ^ Hu{x)H^{y), (2.1) 

k=0 
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2.2 Motivation of Conjecture 3 


2 Theorem 2 and Conjecture 3 


built out of normalized degree k Hermite polynomials Hk, orthonormal with respect to the weight e ^. Define 

Ga^Ak) by 


If ” ' 

ga^m )=- K?'))=Y} n n' 

” „ „ l^iKj^n i=l V 


e X 


1 , 


Xj < do 


1 - Xj > do 


X 6.x j 


By (1.35), we have 

with = I - . 

Similarly as in (1.41), we have 


Gio,„(/c) = M^„,„(log(l -k)) 




//„(4o,0) 


( 2 . 2 ) 


(2.3) 


Ga,M) = M,„,„(log(l - ^2)) = _£(1 - K^fEAKAo) = det(l - 


(2.4) 


k=0 


where ^Ai +,, 0 ) is the integral operator with kernel Kn acting on [do, + 00 ), and the determinant is the Fredholm 
determinant (the proof of the last equality in a more general setting can be found in [28, §23.3]). 

Another well-known result is the convergence of the kernel K„ to the Airy kernel 

Ai (x)Ai'(y)-Ai(y)Ai'(.v) 


Kpa{.x,y) 

if ;c, y are scaled properly around V^: 


x-y 


(2.5) 


:- kA -I- — -, -I —— -I = Kpa{u,v) + e 








-c(|m|+|v|) 


0 ( 1 ), 


( 2 . 6 ) 


uniformly for u,v > -M, M > 0, for some c > 0. Using a slightly stronger version of this Airy kernel limit, as 
in [16], one shows the convergence of the associated Fredholm determinants: if we scale do as in (1.7), we have 




(2.7) 


uniformly for t e (-M, -too) for any M > 0, where ^Ai is the integral operator with kernel Kai acting on 
L3{t, -too). 

Using the Tracy-Widom formula (1.15) together with (2.3) and (2.7), we obtain 


(do, 0 ) exp ^ (r - t)u{T-,l3f6)^ (1 + o(l)). 


as R ^ 00 . This proves (1.13). 


( 2 . 8 ) 


2.2 Motivation of Conjecture 3 

In the case where do = d with d 6 (-1,1), asymptotics for the Hankel determinants //„(do,/3) were obtained 
in [21] and are given by (1.5). The dependence of the error term on d was not made explicit in [21], but it can 
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3 RH analysis of orthogonal polynomials 


be seen from their analysis that the error term in (1.5) gets worse if A approaches +1. We hope that by a careful 
inspection of the estimates in [21], one can strengthen the error term and obtain 

H„(A V^,y8) = (do, 0) G(1 + I3)G{\ -13){\ - x 

X exp(2m^(arcsind + dVl -d^))|l |Re^| < (2.9) 

for some y > 0. The error term must be uniform as d T 1 at a sufficiently slow rate such that n^^^(l - d) is 
sufficiently large, say larger than some fixed M > 0. 

We now take d = 1 + jl with -t > 2M. On the one hand, we can apply (2.9). Expanding the right-hand 
side of (2.9) for large n, we obtain, after a straightforward calculation, 

log Hn(Ao,l3) - log Hn (do, 0) - ninjA 

= - \0'log{-t) + log(G(l -P)G{1 +[A)) - 3[A^\og2 + e„{t\ (2.10) 

where \e„{t)\ < cl\tY + d for some c,d,y > 0, if n and -t are sufficiently large. We thank the referees 

for pointing out the n-dependence of this error term. 

On the other hand, by (1.13), we have 

log//„(do,/3) - logid„ (do,0) - ninfi = logdet(l - |+ o(l), n ^ oo. (2.11) 

Comparing (2.10) with (2.11), we obtain 
logdet(l-/c2i^Ai|[,„^)) 

= -^i/3(-tf^-^/3Hog{-t) + \ogiG(l-/3)G(l + fi)) - 3/3^ \og2 + en(t) + o(lX (2.12) 

as n ^ oo. Letting first n ^ oo and then t -oo, we obtain (1.17). The total integral identity (1.18) now 
follows easily from (1.17) and (1.15). 

3 RH analysis of orthogonal polynomials 

3.1 Overview of transformations 

Following [19] (see also [11] and [24]), consider the RH problem for the matrix-valued function Y{z) analytic 
in both upper and lower open half-planes with the following jump condition on the real axis: 

y4^) = y_(x)|Q (3.1) 

where y±(x) is the limit of Y(x) as z approaches x from the upper (-I-) or lower (-) half plane, and with w{x) 
given by (1.2). Y has the asymptotic condition 

Y(z) = |/ + as z ^ cx), (3.2) 
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3.1 Overview of transformations 


3 RH analysis of orthogonal polynomials 


where cr 3 is the third Pauli matrix 


0-3 



The explieit solution of this problem is 




Pniz) 

-2mh-\pn-Y{z) 


{Ini) 
-h 


-1 f Pn(x)w(x) ^ 


_ -00 - z 

1 p„-i(x)wix) 


n-l 


dx 

dx 


(3.3) 


where p„ and are the monie orthogonal polynomials of degree n and n-l with respeet to the weight 
w(x) = w(x',Ao,/3) defined in (1.2), and hn-\ = /_^^p„-i(x)^w(x)dx. 

This RH problem for Y has been studied asymptotieally, for large n and with To sealed as in (1.7), in [33]. We 
give an overview of the series of transformations eonstrueted in this asymptotie analysis, but refer the reader to 
[33] for more details. Define the function T(z) as 

r(z) = z6C\R, (3.4) 

where ^ 

/ = -1 - 2log2, g(z) = [ log(z - 5)(/f(5)d5, z 6 C \ (- 00 ,1], ilf(s) = - Vl - 52 . (3.5) 

7-1 n 

Here, the logarithm is in its principle branch with branch cut in the negative direction, and i^{s) > 0 on (-1,1). 
As usual, this g(z) satisfies certain variational relations: 

g+(z) + g_(z) = 2z^ + I, z 6 (-1,1), (3.6) 

g+(z) + g-(z) < 2z^ + l, z6]R\[-1,1]. (3.7) 

Additionally, its jump across the real line is described by 

f27r/, z < -1, 

,?+(z) - g-(z) = 1 2m d^, z 6 [-1,1], (3.8) 

[ 0 , z>l. 


Let if/{z) be the analytic continuation of ij/ into C \ ((- 00 , -1] U [1, 00 )). Introduce the function h{z) as follows: 


h{z) = -ni j i/f(y)dy, z 6 C \ ((-00,-1] u [l,oo)), 
and define a piecewise analytic function S in lens-shaped regions (see Fig. 1) as follows: 


(3.9) 


5(z) = T{z) ■ 


/, 

^ 1 0 
—g-OT/3g-2n/i(z) j 

1 0 ^ 


outside the lenses, 

, in the upper half-lens, 

in the lower half-lens. 


(3.10) 


As shown in [33], the function S has jumps on the lens-shaped contour shown in Fig. 1. As n ^ 00 , the jump 
matrices tend to the identity matrix everywhere except on (-1,1) and in small disks U~^ and around -1 
and 1. To obtain asymptotics for 5, an outer parametrix and local parametrices near -1 and -l-l have to be 
constructed. 
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3.2 Outer parametrix 


3 RH analysis of orthogonal polynomials 



Fig. 1: Opening of lens, case f < 0. 


3.2 Outer parametrix 

For z outside small disks around -1 and +1, S can be approximated for large n by an outer parametrix 
which is analytic except on [-1,1], tends to the identity as z ^ oo, and has the jump relation 



(3.11) 


It is given explicitly (see e.g. [21]) as 



(3.12) 


where 



(3.13) 


The branch of uq is chosen so that ao(z) ^ 1 as z ^ oo. 

3.3 Local parametrix near 1 

In order to obtain asymptotics for S also in neighborhoods of +1, local parametrices have been constructed in 
[33]. Near -1, this local parametrix was built using the Airy function, but we do not need its explicit form. 
Near +1, it was built using a model RH problem associated to the Painleve XXXIV equation. 

The local parametrix P^^\z) is analytic in t/', except for z on the jump contour for 5, and it has the same 
jump relations as S for z on the jump contour for 5, inside t/'. On the boundary dU^, it satisfies the matching 
condition 

P^^\z) ■ P^°°\z)~^ = I + asn ^ oo, uniformly for z 6 dU^. (3.14) 

It takes the form 

F(‘)(z) = F(z)0(4'(z); 

where E is an analytic function in t/\ d) will be specified below, and ^(z) is a conformal map near 1. The 
conformal map ^'(z) and the parameter r are given by 



(3.16) 


Here, the multivalued power is the principle branch in C \ (-oo, 0). The Taylor expansion at 1 is 



(3.17) 


(3.18) 
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3.3 Local parametrix near 1 
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The analytic pre-factor E can be expressed as 


and 0(^; r) is given by 

W;t) = (q yjo(4' + T;T), (3.20) 

where 'To(^; t) is the solution to the following RH problem. 



Fig. 2: The RH problem for ToC^). The rays meet at ^ = 0. The union of the 
rays is referred to as Fiy^. 


'To is analytic off the contour shown in Fig. 2 and satisfies the following jump and asymptotic conditions: 


'To,(^) = 'To_(^) • 


0 1 
^1 0^ 

vl h 
‘ 0 1^ 
-1 0 


/ 


^6 71 , 

^ 6 72 U 74, 

^6 73 - 


(3.21) 


'Po(^) 

'Po(^) 


m 

a b 
c d 




oo. 


(3.22) 

(3.23) 


where the matrix elements of m as well as a, b, c, d are some functions of r, k is given by (1 . 16), and M is a 
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3.4 Lax pair for 'Fq and the Painleve XXXIV equation 


3 RH analysis of orthogonal polynomials 


piecewise constant function defined as follows 


^1 0^ 

vO h 
^ 1 0^ 




V 1 1 

0 1 


^6/, 

^6/7, 

^ 6 HI, 
^elV. 


(3.24) 


All multivalued functions above are in their principle branches with branch cuts along the negative half axis.'Po 
is uniquely determined by the above conditions. Note that all higher order terms in the expansions of 'Pq 
also functions of r. 

The function defined in (3.15) is the same as the one in [33], but it has to be noted that our functions ¥0 is 
defined in a slightly different way compared to [33], which will be convenient later on. We have the relation 




2-^3/6»j/o = 2^/^; T = 5). 


(3.25) 


where denotes the solution to the model RH problem of [33]. 

By (3.20) and (3.22), it is straightforward to verify that O admits the asymptotic expansion 


m 


d>{C,r) = \I+ — +0 


1 




.^3/4 1 / 1 -'h-k’®.43 


VaH 1 

where we have the following relation between m = m(r) and m® = 


as^ 


m = 


<i> , ^ ^ <i> 

mn + ^rn,, + --- mi2-- 

< 1 . ^ 

^21 + -J 


mfi + 


16 


21 


ir^ ir^ ^ 

:m,, - 

a, a, 

"*22 - ^"*21 - Y + ^ 


(3.26) 


(3.27) 


3.4 Lax pair for To and the Painleve XXXIV equation 


From the RH conditions for 'Pq: there is a standard procedure to deduce differential equations with respect to 
the variable ^ and the parameter r. Here, our approach deviates from the one in [33]. 

Consider the functions U := ^'Pq^ and V = ^'Pq'. Because the jump matrices for 'Pq are independent of ^ 
and T, U and V are meromorphic functions of Using the behavior of'P q at infinity and 0 given in (3.22) and 
(3.23), we obtain after a straightforward calculation that 'Pq satisfies the Lax pair 


S^o 

dr 


(^;r) 

(e,r) 




mT) = V(^;T) + 





-ir/2\ 

2;rz^ 

0 1 

1-17121 

2mn\ 

1-1 

m2i / 



(3.28) 

(3.29) 
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3.4 Lax pair for 'Fq and the Painleve XXXIV equation 


3 RH analysis of orthogonal polynomials 


where a, b, c, d and the matrix m, which are functions of the parameter r (and also of P), were defined in 
(3.22)-(3.23). We can also refine the expansion for as ^ ^ oo: 




Om-l 


dr 


^-P-V{PT) = 


1 dm i 
^ dr ^ 


m. 


0 0 

-1 0 


m. 


0 1 
0 0 


0 1 
0 0 


m,m 


+ 0 


(3.30) 


Since this expression obviously has to be zero, equating its (21) entry to zero gives us the useful relation 


1 2 * / 
"ill = 2"*21 - 2"*21- 


(3.31) 


Note that m^j is the square of the matrix element m 2 i. The compatibility condition of the Lax system (3.28)- 
(3.29), 


V,-Ur = [U,V], 


(3.32) 


becomes 

(0 1 

0 0 


+ I 


-"*21 2 m ; 1 

0 m( 


21 




- (ac)' 

-(c^) (ac)' 


rm 2 i\ . (c^ -2ac\ ik^ l2pmi\-a^ 2 a^m 2 i - 4acm 

2ni^ \ 2 Pm 21 - 2ac cd - 2pm\i 


2ni \ 0 


(3.33) 


This equation can be separated into two equations for each power of p From the resulting system one can 
extract the equations 


« 2 w ^ 




hr-'Pi = K)' 


and 


(1 -I- 2iTm2\ — 4m; j)^ -I- 4 (2m;j - ir) {2im"^ -I- 2/mn (r -I- 2/m;j) -I- rm^j -I- m2\) = 0, 
which, with the help of (3.31), reduces to 

1 -I- 32T{rn2if + 32i{m2if + dzm^'j - 4(m2j)^ - Airm'P + Sm^j (m'P - ir^^ = 0. 
This equation is a disguised version of the 34th Painleve equation for the function 

y(T) = -zm;i(T) - ^ = -/(m®) (r), 

namely. 


4y^ -I- 2Ty + 


(yrY 

2y 


Equation (3.34) also provides us with another representation of y(T): 


y(T) = zTim 
f-»o 


d'Fo(^) , 

(^) 


(3.34) 

(3.35) 

(3.36) 

(3.37) 

(3.38) 

(3.39) 

(3.40) 


J21 


Moreover, from (3.31), we obtain an additional expression for y which does not involve derivatives: 


y(T) = 2mii(T) - mj^ir) - t/2. 


(3.41) 
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3.5 Proof of Theorem 1 


4 Asymptotics of m(t; k) as t ^ +oo 


3.5 Proof of Theorem 1 

In [33, Corollary 1], it was proved using vanishing lemma teehniques that the RH problem for s) is 

solvable for all real values of 5 if jS is sueh that |arg < n, and thus for all p sueh that | RejS| < 1 /2. Beeause 
of the explieit relation (3.25), this implies that the RH problem for 'Pq is also solvable for all real values of 
T if I Re/3| < 1 /2. This in turn implies that the funetion y(T) = y(T; jl) defined in terms of 'Pq by (3.41) is 
well-defined and eannot have singularities for real r if | Re/3| < 1 /2. 

If we define m(t; k) by m(t; k)^ = y{T;/3) with = I - then it is easily verified by (3.39) that u solves 
the Painleve II equation (1.8). By exploring the asymptotie behavior of y(T;/3) (or, equivalently, of u(t;k)) 
as T ^ + 00 , we will be able to identify u{t; k) as the Ablowitz-Segur solution of the Painleve II equation 
eharaeterized by (1.9) and (1.10). This identifieation, whieh will follow from (4.33) below, eompletes the proof 
of Theorem 1. 

3.6 Final transformation 

Introduee the new funetion 


R(z) = 5(z) • 




Z6C\t/-lUt/lU(-l,l), 

Z 6 U^-^\ 

Z 6 


(3.42) 


which tends to the identity matrix as z ^ oo and whieh has jump matriees Gr on the eontour T« that tend to 
identity as n ^ oo: 


G«(z) = 


y~)^z)(p[“>(z))“', 

. p(oo)(z)|p(-l)(z))-\ 

p(oo)(z)|p(l)(z))-\ 


Z6 (-1,1)\ t/-l U U\ 

zedU-\ 

zedU\ 


I + O 


1 


-I- \z\)P 


z 6 T^ 


(3.43) 


This, in turn, implies (see [13]) that for suffieiently large n 


R(z) = 1 + 0 


’ uniformly for z 6 C \ Yr, 


where T^ is the jump eontour for R. 


(3.44) 


4 Asymptotics of u(t; a:) as t ^ +oo 

From Seetion 3.4, we know that y(T;fi) defined by (3.38) solves the Painleve XXXIV equation (3.39), and 
this implies that u defined by u{t; k)^ = y(T; /3) (with relation (1.16) between k and /3) solves the Painleve II 
equation (1.8). We now proeeed with proving the asymptoties of y(T; yS) = u(t; k)^ as stated in (1.9). In this 
seetion it is assumed that t > 0. The analysis performed here is largely analogous to the one in [22] with one 
additional new teehnieal feature - the need to introduee an additional triangular parametrix near the point z = 0 
of the diseontinuity of the (triangular) jump matrix (see Subseetion 4.3.2). 
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4 .1 Rescaling and shift of the jump contour 


4 Asymptotics of m(t; k) as t ^ +oo 


4.1 Rescaling and shift of the jump contour 

Introduce 

A(z) = t). 

One can easily cheek that it satisfies the following RH problem. 

(a) A : C \ Fip;, ^ is analytie. 

(b) A has the same jump relations as 'Pq- 


(e) A{z) = (l + 0(^): 


,0-3/4 j_ 


1 -i 


pi2Qp/2^^l/2y^ 


V2 1 

(d) A(z) inherits its behavior at z = 0 from 'Pq very easily. 
From (3.40) we get 


as z ^ 00. 


y(T) = —— lim 


-\/t [ dz 


,-i 

2 —;— Aiz) 


J21 


We shall further write 




With respeet to the domains defined in Fig. 3, define 


B(z) 


Aiz) 


Aiz)- 

Aiz), 


1 0 

vl h 
‘ 1 0^ 
-1 1 


z 6 ir, 

, zeiir, 

z 6 7 u ir u nr u iv. 



Fig. 3: The contours F b and the RH problem for B{z). 


This funetion satisfies the following RH problem. 

(a) B : C\Yb ^ is analytie. 
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(4.1) 


(4.2) 

(4.3) 


(4.4) 












4.2 Normalization at infinity 


4 Asymptotics of m(t; k) as t ^ +oo 


(b) B^(z) = B_(z) 


1 


,0 1 , 
^1 0^ 


vl h 
^ 0 1^ 
v-1 Oy 


0 1 


Z 6 TbI, 
z 6 762 U 764 , 
Z 6 7B3, 

Z 6 7B5- 


1 -i 


-<| z 3/2 


(c) Bfe) = (/ + 0(l))z"«^^_. , 

(d) B{z) has logarithmic behavior near z = 0. Namely, 

sfe) = gfe)(/ + A(o ; 

where B{z) is some analytic function, M+ = I and M_ 
branch with branch cut along the negative half axis. 


as z ^ oo. 


logz M+, Z6 


(4.5) 


0 


,-2nip 

1 


The logarithm is in its principle 


Finally, the expression for y{T) remains unchanged compared to (4.2), 

mz) 


y{T) = 


— lim ^ . 

yr z^o [ dz 


-B{z)- 


(4.6) 


J21 


This transformation is an important precursor of the introduction of a new g-function that will allow us to 
“undress” the behavior of the RH problem at infinity (for similar transitions see e.g. [22] and [32]). 


4.2 Normalization at infinity 

We now introduce the following g-function, 

g(z) = ^(z + l)^^^ -n < arg(z + 1) cr. (4.7) 

Note that 

g(z)-(^z'^' + z'^") = asz^oo, 2 ^ (- 00 ,- 1 ]. (4.8) 

Next, define 

C(z) = |q (4.9) 

The constant prefactor in this definition is needed to conserve the <9 (^) term in the asymptotics as z ^ oo. C 
satisfies the following RH problem. 

(a) C : C \ r6 ^ is analytic. 
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4.3 Construction of parametrices 


4 Asymptotics of m(t; k) as t ^ +oo 


(b) CAz) = C_(z) 


-1 0 

J g-2ig(z)' 

0 1 






0 

1 0 

«2ig(z) 1 

^ 0 1^ 


Z 6 TbI, 
z 6 752 U 754, 

Z 6 753, 

Z 6 755. 

as z ^ OO, z ^ r5. 


(d) 


C(z) = C(z) ^ |q qI log zj for z 6 C- near 0, where C is analytie in a neighborhood of 


0 and M+ are the same as before. The braneh of the logarithm is ehosen as before. 
From the definition of C, 

' d5(z) 


dC(z)^. 3-1 
z— 5 —C(z) 
dz 


J21 


dz 


-B(zy 


+ 


J21 


sg(z)[zB(z)o- 3 B(z) ^ 
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The seeond term in this expression tends to zero as z ^ 0 due to the behavior of B{z), henee 


7(t) = —— lim 


■\/t 2^0 [ dz 


dC(z)^. 3-1 
— C(z) 


(4.10) 


(4.11) 


J21 


4.3 Construction of parametrices 

4.3.1 Global Airy solution C*^' ^ 

The jumps of C(z) near z = -1 are very similar to the jumps of the standard Airy RH problem. Let us look for 
a funetion that satisfies the following RH problem. 

(a) : C \ Ts ^ is analytie. 

/1 e“2ig(z) 

(b) C'^‘yz) has the same jumps on r 5 \ [-1, +oo) as C(z) and its jump on (-1, +oo) is L ^ 


(c) C(A‘>(z) = (/ + 0(j))z‘^^^^^ I as z ^ cx,. 

We seek in the form 

C(Ai)(z) = C(Ai)(2)5^tfeV3_ 

If we define an auxiliary matrix funetion (whose jumps are shown in Fig. 4) 

yo -y2 


yo yi 

[yo y'lj 


(4.12) 


' 

( \ 
-yi -yi 

in II, 

/ 

= 


ry'i -y'l. 


v 


( \ 


< 

/ 


-yi yi 

in III, 



-3^2 3^1. 



V 




in I, 
in /y. 


(4.13) 
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4.3 Construction of parametrices 


4 Asymptotics of m(t; k) as t ^ +oo 


where 

= Ai(z), yi(z) = y2(z) = 

then a standard argument shows that must have the form 

C(^')(z)= V^l^. 1)), 

This in partieular implies a refined asymptoties for 




m 


Ai 


+ 0 


.0-3/4 


1/1 -i 


V2H 1 


m 


Ai 


cr3 li to 1 
0 


Airy solutions like this are diseussed in mueh detail, for example, in [14]. 



Fig. 4: The standard Airy RH problem. The union of contours is referred to as 
Fa.. 


(4.14) 


(4.15) 


(4.16) 

(4.17) 


4.3.2 Local solution 

We also need a loeal parametrix for C near z = 0. Let [7° be a small open disk around 0 of radius less than 1, 
say, 1/2. Then we have to find the funetion which satisfies the following RH problem. 

(a) : t/° \ [ 0 , + 00 ) ^ is analytie. 

(b) C+ (z) = (z) Iq ^ I, z 6 ( 0 , + 00 ) n [7® (eontour oriented to the right). 

(e) C*^®^(z) = / + o ( 1 ) as 5 ^ 00 , uniformly on dU^. 

(d) C'^®^(z) ~ ^ |q qI logzj as z ^ 0. The braneh eut of the logarithm is along the positive half 

axis. 
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4.4 Final transformation 


4 Asymptotics of m(t; k) as t ^ +oo 


Due to the simple algebraic structure of the jumps, this problem ean be solved exactly in terms of integrals of 
elementary functions. Namely, the solution is 






C^°^(z) = 

2;rz /o 

^-dz' 

z - z 

, Z6C\ 

0.^ 


lo 

1 J 




(4.18) 


This function clearly has the requested jumps and has the same general logarithmic behavior near z = 0. 
Moreover, this funetion satisfies the matehing eondition on dU^ and, in fact, with some c > 0 we have 

= I + 0{e~‘^^) as 5 ^ oo, uniformly on dU^. 


We will also need the faet that 


limz 




dz 




j = 

2ni 


o 

o 


(4.19) 


(4.20) 


D{z) 


z 6 C \ t/o u Tfi. 


(4.21) 


4.4 Final transformation 

Using the functions built in the previous subsection, we can now perform the final transformation of the RH 
analysis in the case where r ^ +oo. 

Define 

fc(z) ■ (C(«>(z))“' ■ (C(^‘^(z))“', z 6 t/o \ R, 

[C(z) ■ (c(A‘)(z))“', 

This function has the following properties. 

(a) D : C \ +ooj u ^ is analytie. 

(b) Assuming the countercloekwise orientation of dU^, 

[o 1 , 

[c(A‘)(z) ■ (C(0)(z))“' • (c(A'^(z))“', 


D^z) = D_(z) 


Ci^‘\z) 


■(Ci^'^(z)) \ Z6(5,+oo), 

Z 6 5t/°. 


(4.22) 


(e) D(z) = I + as z ^ oo. 

Using the asymptotie expansion for it is easy to eheek that, with some c > 0, 


as 5 ^ oo, uniformly for z 6 ( -, +oo |. 


1 


(4.23) 


(D_(z))“'D+(z) = / + <9(e-"l^l) 

As for the jump on dU^, by virtue of (4.19) and the boundedness of it is also elose to the identity matrix: 

(D_(z))-' D+(z) = C(^')(z) • (I + 0{e-^^)) ■ (c^^'^z))"' = / + (9(e-"0 (4.24) 

as 5 ^ oo, uniformly for z & dU'^, with some c > 0. 

Using these estimates, in a standard way [13] one shows that, for any z 6 C \ To, 

e" 


D(z) = 1 + 0 


1 + kl 


as 5 ^ oo, c > 0. 


(4.25) 


The error term is uniform on eompaet subsets of C \ T/). 
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4.5 Asymptotics fory and uniformity of error terms 


5 Asymptotics of the recurrence coefficients 


4.5 Asymptotics for y and uniformity of error terms 


Following the transformations (eqs. (3.20), (4.1), (4.4), (4.9), (4.21)) 

backwards, we can recover the connection between the asymptotic expansions of O and D. Namely, for large 
z 6 C \ (//' U III'), we have 


0(2) = - 1) 


0 -1 
-i 0 


r'^3/4o(Ai)(z). 


Next, we write, as usual. 


D(z) = /+ —+ ^ 

z z^ 


+ <9 


as z ^ oo, 




0(z) = 7 + — + 


+ 0 


1 


, 0 - 3/4 


1/1 -i\ _2 


V2H 1 


e ^ ^ as z ^ 00 , 


(4.26) 

(4.27) 

(4.28) 


z-^ 

Let us refer to the matrix coefficients in front of z~^ in these expansions as and (they are only functions 
of r). Using (4.26), we see that each of the matrices in the expansion for O is merely a linear combination 

of a finite number of the matrix coefficients m^’*, with coefficients rational in 
Using (4.25), one shows that the matrices are exponentially small, 

DJ 


m‘^’*(r) = 0(^e as r ^ + 00 , with some c > 0 
for all k. It immediately follows that 

m^’^(T) = as r ^ + 00 , 


(4.29) 

(4.30) 


thus m®’* are bounded at large r. 

These facts imply that the asymptotic expansion (3.26) for O is uniform for t 6 [tq, + 00 ) for any tq 6 R. 
Since C(z) = D{z) ■ C*^aO(^) . c*^°^(z) in 77° and both D{z) and are bounded there, we get from (4.11) that 


yfr) = —— lim 


z^o [ dz 


clC(z)^, ,_i 
—C(z) 


J 21 


;f2g-2ig(0) 

In yfr 


D(0)C(Ai)(0) |0 1| (^(0)C(Ai)(0))“' 


for large r > 0. 


J 21 


C^Ai)(0) = 




+ 00 . 


(4.31) 

(4.32) 


From (4.12), (4.15), and the asymptotics for O'^ao it follows that 

0 - J_/l -M 

0 / V 2 \-* 1 

Taking into account that D(0) converges to I very rapidly and using the definitions (4.7) and (4.3), we arrive at 
the final expression for the asymptotic behavior of y(T;l3), where we emphasize the dependence on fS: 


yir-,/3) 


43/2 / K 

" -2 * I 




Together with y(T; /3) = u{t; nf, this proves (1.9). 


5 Asymptotics of the recurrence coefficients 

In this section, we compute asymptotics for the recurrence coefficients 7?„ and Our calculations in this 
section are similar to those in [33], but we believe it is convenient for the reader to give some details of the 
calculations because of differences in notations. 
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5 .1 Auxiliary asymptotics of Gr 


5 Asymptotics of the recurrence coefficients 


5.1 Auxiliary asymptotics of Gr 

We now need to compute the precise asymptotic behavior of Gr, the jump matrix for R (see (3.43)). Finding an 
explicit expression for the two leading terms in Gr on dU^ is the most sophisticated part of this calculation. 
First, expand 


Gr{z) = P^°°\z){P^'\z))' = EizT^ forz 6 dU^ 


(5.1) 


with E{z) defined in (3.19). Recall that we have the asymptotic expansion (3.26) for O, uniformly for r > tq 
with any fixed tq 6 R. Therefore, one verifies using (3.12) that, as n ^ oo. 


Gs(z) = e“'’"'^2_(' ‘WJ—L\ 


■n\> mz+il 


0 - 3/4 


I- 


m 
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X 


z-ir”" 1 /I -A - 


iTT/jcrj/l 


\Z+1I 1 

which gives us the following expansion of Gr as n ^ 00 : 

Gr(z) = / - Gi{z)n~^'^ + G 2 iz)nr^^^ +0{n~''^ , uniformly for z 6 dU^ and r > tq, 

where 


Gi(z) 


zm?, (z + 1)'^^ 


, 1/3 


iz-\f'\{zyi^ 


1 


-le 


-inf) 


-1 


and 


G2{z) = 


I 0 e™-® 


^(z) 0 /- 


(5.2) 

(5.3) 

(5.4) 

(5.5) 


5.2 Asymptotics of R„ 

We can use the following simple identity for the recurrence coefficient defined in (1.4) (see e.g. [11]): 


Rn = zn[2ZW2i. 


where the matrix is defined in terms of the large z expansion of Y: 


Y(z) = \I + 


+ (9|liu«^3 


rrf'it) 


In order to compute zn^, we will need similar large z expansions for the following functions 


R(z) = 1 + 

p(‘zo) — I ^ 


z 

m^(t) 


+o\G. 


-t-oK . 




(5.6) 


(5.7) 


(5.8) 

(5.9) 
(5.10) 
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5.2 Asymptotics of R, 


5 Asymptotics of the recurrence coefficients 


Unfolding the transformations Y T S i-^i?at large z, we obtain the identity 


( 2 „)«^ 3/2 (rn^ + ( 2 r )-”^ 3/2 ^ 




Since we can reformulate the RH problem for R in terms of an integral equation 


R_{Z) = I+ — 
2ni 


1 f R_(z)(G«(z')-/), , 


'Tr 


z -z 


dz', 


we have 


= R_iz'){GR(z')-I)dz'. 

2m Jr, 


(5.11) 


(5.12) 


(5.13) 


Next we iterate the integral equation to find an asymptotic expansion for R_ as n ^ oo. Given that integration 
over the contours other than dU^ gives only a O contribution (because the jump matrix Gr = I + 0{n~^) 


on F/; \ dU^), the large n expansion for is 

^(Gr(z ) -/)dz' - ^ ^pi{z){GR{z)-I)dz +0{n~^). 

00 00 

We can now substitute the asymptotics (5.3) to get, after a straightforward calculation, 

res Gi(z) - res G 2 (z) + res Gi(z) • res +0(n~^). 
as n ^ oo. Now, from (5.4) and the expansions (3.17), (3.18) we find that, as n ^ oo. 


(5.14) 


(5.15) 


resGi(z) 

Z=l 


and 


as well as, from (5.5), 


res 


:=1 Z - 1 


res G 2 (z) 

Z=l 


zm® (0 ( 

' 1 

-ze“^\ 

2 1 


-1 ) 

zm® (0 / 
10 \ 

1 

-1 / 

( t ) 

1 ^ 


2 


or 




(5.16) 


(5.17) 


(5.18) 


Note that res 2 =i Gi is nilpotent, thus the third term in (5.15) is negligible. Furthermore, from the relations (3.27) 
and (3.31), we find 

1 / rt, \2 i 


<i) ^ t V 

- 9 Kij 


Substituting all previous results into (5.15) we obtain the final formula 


m 


^ _ im^i(0 I 1 0 e 




rjV 

2 


inf} 

e 0 


+ 


0(n'). 


(5.19) 


(5.20) 


The second matrix in (5.11), m“, can be easily found from (3.12): 


i I 0 


m = 




2\-e-'^P 0 /■ 


(5.21) 
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5.3 Asymptotics of Pn (/lo) and Q, 


5 Asymptotics of the recurrence coefficients 


All the operations performed to obtain the asymptotics of from the asymptotics of O preserve the uniformity 
in T 6 [tq, oo), or equivalently t e [to, oo), for any tq, to ^ R. 

By substituting and m“ into (5. 1 1) we find the large n expansion for R„, 


Rn = 


2 n (^21 + ^21 ) 

/m® ■ " 

\ 2 2 


_ 2/3 I 

—-n + 




tfl^ ini^ 7 


+ 0(n-') , (5.22) 


which, by (5.19), simplifies to 


Rn{Ao;/3) = ^ + ( 9 ( 1 ) as n ^ oo, for all t 6 R, uniformly for t > to. 


(5.23) 


since = y. This result holds for all t 6 R and jS such that | Re/3| < 1 /2. This asymptotic series formally 

matches the classical Hermite recurrence coefficient asymptotics when T > 1 (t ^ +oo) and the non-critical 
asymptotics from [21] when T < 1 (t ^ -oo). 


5.3 Plancherel-Rotach type formula and asymptotics for Qn 

We can express the orthogonal polynomial p„ in terms of the RH solution Y, 

Pn (t V^) = Tii (d V^) = lim {2nf'^ 5 „ (z) 


(5.24) 


where the limit for S is taken for z approaching 3=1 + from the upper half plane, outside the lens-shaped 
region in Fig. 1 . If z lies in this region and z & (the small disk around 1 in which the local parametrix 
was constructed), then we can unwind the transformations S O 'To to obtain 


S(z)=R(z)P^'Hz)=R(z)E(z) 


1 


IT 




0 1 


'To (z) - t) e' 


-m/jiT3/2^lPzf^^cr3 


(5.25) 


In order to compute the limit z ^ 3, we need to use the small ^ expansion of 'Tq in sector I given in (3.23). 
After a straightforward calculation, using also (3.44) and (3.18), we get 


Z—*/l ' \ // 


Therefore, 


Pn (a V^) = (2nf^ (l + O . 

From the Lax pair identity (3.34) and = u{f, k)^, 

2nu (t; a:)^ 


OO. 


OO. 


and 


{ic{np)y 


= ± 


(t; k) 


(5.26) 

(5.27) 

(5.28) 

(5.29) 


The right hand side does not depend on the sign of k (indeed, changing k to -k changes u to -u), and we can 
verify which sign is correct using the asymptotics for c as t ^ oo. Since 


c(t) = fim('To(z))2i, 

z -»0 
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5.3 Asymptotics of Pn (/lo) and Q, 


5 Asymptotics of the recurrence coefficients 


working backwards along the transformations A B C D for both r ^ +oo, we can easily 

recover the asymptotics for c (r). It turns out that 

_ 2 ^ 3/2 

ic (t; P) = - (l +0 as t ^ +oo, (5.30) 

V2ri/4 

which implies that the correct sign in (5.29) is +. 

Remark 12. In the special case = 0, the model RH problem for ¥0 reduces to the Airy model RH problem. In 
this case, we have 

ic{T;f3 = 0)= V^Ai(T), (5.31) 

which is indeed consistent with (5.30). 

We thus have 

n ^ 00. (5.32) 

K 

Now we need the expansion of g{z) near z = 1: 

g(z) = ^ - log2 + 2(z - 1) - ^2^/2 (z - + 0(z - I)' as z ^ I. (5.33) 

Substituting 3=1 + we have 

4 

2ng+ (A) = n - In log 2 + + -is + O as n ^ 00. (5.34) 


This gives us the asymptotics as n ^ 00 of the polynomials near the critical point, 

(3 . (5.35) 

By multiplying the recurrence relation (1.4) by pn{x)w{x) and integrating, we find 

Qn =-hTn pniyA^j2^ sivA\{inp). (5.36) 

Note that 

hn = - lim 2m ¥21 (z V^) (z V^) = - lim 2m (2n)” e"^z S i 2 (z) = 

= -2m {2nT e"' (m“ + m’l^ , (5.37) 
thus the following large n asymptotics hold for the normalizing coefficients 


nyflnn 


= on n “ /mfi(0n - mfi(0n +<9(n ^)) 


2"e 

This proves (1.26). Equivalently, this result can be deduced from the identity 


hn — 




(5.38) 


(5.39) 
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6 Hankel determinants: alternative proof of Theorem 2 


expressing /?„ as a ratio of two Hankel determinants, together with the asymptoties (1.13). Substituting (5.38) 
and (1.27) in (5.36), we obtain (1.25). 

Lastly, we note that we ean easily obtain the asymptoties of the eoefficients in (5.38) as t ^ -oo. Formally 
this ean be done by computing an antiderivative of the asymptotics of y (t). The following asymptotics were 
obtained rigorously by solving the RH problem in the limit t -oo. For |RejS| < 1 /2, 


im. 


21 


(t;jS) = 


1 ( ni 

4i(-T)\ Tifi) 


r(l+y8) . 
r(->s) "" 

30^ 


2(-t) 


-mr-p) 


+ 0 


1 


(-T) 


5/2-3|Re,8| 


as T ^ -oo. (5.40) 


When yS = Ik, /<• 6 R, this becomes 


I— K / 4 

im® (t; ik) = 2k h- 


2(-r) + 3/dog(-T) + 6/dog 2 - 2 arg F(//c) I + 

3/c^ / 1 

+ ;r-—- +0\ -^1, asT 

2(-t) \(-t)^^^ 


-oo. (5.41) 


When 13= 1/2 + iy, y e R, we have 

1 


/m® |t; 2 + /yj = V^|2y-tg|^jJ+(9p|, as r 


(5.42) 


These formulas complement Theorem 5. 


6 Hankel determinants: alternative proof of Theorem 2 

6.1 Differential identity 

Here, we derive a differential identity for the logarithm of the Hankel determinant Hn(Ao,l3). It is expressed in 
terms of Y defined in (3.3). 

Proposition 13. We have 

— logH„{Ao,/3) = - sinnfi (6.1) 

oAq n ^ ''21 

Here ' is the derivative ofY{z) with respect to z- 

Proof. We write Pk = KkPk, Xk = > 0 for the normalized orthogonal polynomials with respect to the weight 

w. We start from the general identity (equation (17) in [27]) 

^ log Hn{Ao, n) = (7^ _ 72), (6.2) 

oAq Kn-\ Kn 
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6.2 Asymptotics for the log. derivative of H, 


6 Hankel determinants: alternative proof of Theorem 2 


where 


Jl= Pn{x)P'„_iix)w{x)dx, 


J2= [ Pn(x)Pn-l(x)w(x)dx. 

Jr 

Here and below dots denote do-derivatives and primes denote jc-derivatives. 
To eompute Ji, we proeeed as follows: by (6.3) and (1.2), we have 


(6.3) 

(6.4) 


■^1 = / Pn{x)P'^_^{x)w{x)dA- [ Pn(x)Pl^_i(x)w(x)dx 

\Jr / Jr 

+ 2ism(n/3)PMP:_Mo)e-^°- (6-5) 

The first two terms vanish by orthogonality, and we obtain 

J: = 2i sm{n/3)(P„P'„_,)(Ao)e-'^l (6.6) 


Similarly, by (6.4) and (1.2), we have 

J 2 = ^\[ P'n{x)Pn-\{x)w{x)dA- [ P'^{x)Pn-\{x)w{x)dx + 2i sm{nlJ)P'^{AQ)P . (6.7) 

OAq \Jr j Jr 

Using the orthogonality relations, we ean eompute the first two terms and we get 

Ji = 2i sin(;rjS) i) (/lo)e“'^“- (6.8) 

Xn-l 

Substituting (6.6) and (6.8) into (6.2), we get 

d ^iK 1 2 

— \ogH„iAo,/3) = -^iPnPn-i-P'nPn-i)(Ao)sm(n/3)e-^o (6.9) 

uAo Kfi 

= -p>„-i)(do)sin(7r/3)e“'^o, (6.10) 

hn-l 

and using (3.3), we obtain (6.1). □ 

6.2 Asymptotics for the logarithmic derivative of Hn 


Let do be of the form (1.7). The results in Section 3 are valid in the limit where n ^ oo, uniformly for t > to 
for any fixed to 6 R. 

Inverting the transformations Y T and T S from Section 3.1, it follows from (6.1) that 

^lnH„(Ao,P) = lim(s-\z)S'(z)) (6.11) 

and the limit is taken in the region outside the lens, see Fig. 1 . By (3.6), (3.8) and (3.9) we have 2g+ (A)-2A^-l = 
2h (A), hence 

^ In H„ iAo,P) = lim (s (z) S '(z)) 


30 








6.3 Expression for r in terms of u 


6 Hankel determinants: alternative proof of Theorem 2 


Near A, we have S = RP^^\ and this implies 

^ In H„{Ao,fi) = ^-smnfi (A) 

dAo ^J7M^ ” ^ ^^21 

+ ^=-sin;r;0 (A)e^'‘^^'^\ (6.12) 

’ h\ 

Sinee R is elose to 7, the seeond term on the right hand side is small. Using the asymptoties (3.44) for R, we 
obtain 

^ \nHM^,l3) = -^-sinnp +Oin-^l^) (6.13) 

dAo ^J^n V ' ^ ' hi ^ ' 

as n ^ oo, uniformly for t > to. To eompute the remaining matrix entry, we ean use (3.15), whieh yields 

Q 11 

— In77„(do,;0) = C{X)—e-^^f^- sinn^ (0) 

oAo V2n ^ ^ 

+ sin;r;0 (o-'(T)£-^(d)£'(d)0(T)) (T) (6.14) 

y/^n ^ /21 V / 

as n ^ oo. By (3.19), the seeond term in the right hand side is of order 0(n~^^^e^h uniformly for t > to. The 
first term will be larger than the last two: by (3.17), we get 

o 1 

— log77„(To,>e) = smnl3 {^o^%) (0) + , (6.15) 

as n ^ oo, uniformly for t > to. Write 

r(T):=(T'-iT';)^^(0;T). (6.16) 

Then, as n ^ oo, 

d 1 

— log H„{Ao,p) = sin (;r;S) r(T) + A {n, t) , (6.17) 

oAo n 

where A (n, t) = uniformly for t > to as n ^ oo. 


6.3 Expression for r in terms of u 


Proposition 14. Let r be defined by (6.16), ¥0 as introduced in Section 3.3, and let u be the Painleve II solution 
characterized by (1.9). The following identity holds. 


d -2m 2 


( 6 . 18 ) 


where k and f3 are related by (1.16). 
Proof. Define 


'Po(^) 


1 -m2i 


'Po(^). 


(6.19) 


\0 1 

This transformation has the advantage that it simplifies the r-equation in the Lax pair (3.28), (3.29). We have 


dr 


'Fol'Po' = 


0 l\ ./ 0 w 

i 


0 0 


-1 or 


( 6 . 20 ) 
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6.4 Proof of Theorem 2 


6 Hankel determinants: alternative proof of Theorem 2 


where w is some unknown funetion of t. In what follows, primes will be used for differentiation w.r.t. r. 
Now, we start from (3.23). In seetor I, we ean write 


for some matrix Ei which is ^-independent. We easily see from (6.16) and (6.19) that 

r(T) = ( 0 ;t) = ( 0 ;t) = £i,2i(t). 

Substituting (6.21) in (6.20), we obtain 

2iW = 

By (5.28), we have 

, 2.7rz 0 

^i2iW = —^M(T;/r)l 

Together with (6.22) and (6.23), this implies (6.18). 


( 6 . 21 ) 


( 6 . 22 ) 

(6.23) 

(6.24) 
□ 


6.4 Proof of Theorem 2 


Asn ^ oo,we haver ~ t,see(3.18). Integrating (6.17) from do = V^(l+to^ ^^^/2)to/li = V^(l+tiR ^^^/2) 
with Iq <ti, we obtain 


ln//„(V^(l + |n“^^^),/3) -ln//„(V^(l + 


h 


■2/3 


),/S) 


- e 


-mfi 


- sin nP 
n 


rh I 

/ r(T)dt -I - —n 

'to V2 


1/6 


A (n, t) dt. (6.25) 


We note that here, as well as in (6.17), 

Writing the left hand side of this expression in an explicit form, one can easily check that there exists a positive 
constant cq such that r > cot for all t > 1 and all n > 1. Hence we can let 0 ^ -i-oo in (6.25) and, taking into 
account that ^ (l + /2^ ,ySj tends to the Gaussian Hankel determinant Hn (do, 0) without the 

jump, arrive at the estimate 


log//„(V^(l + |n-2/3 


1 f°° 

'),/?) - log//„ (do,0) - in/Sn = sinn/S / r(T)dT + , (6.26) 

^ Jto 


or 


//„(V^(i + |r- 


■2/3 


),yS) = V"^”//„(do,0)exp - 




sinn/S 


n 


r(T)dT (1 -to(l)). 


(6.27) 


as n ^ oo. We note that we have replaced dt with dr in the integral r (r) dr. This is justified in the limit 

n ^ oo since r = t(l +0{trT^l^y^ and because r(T), being proportional to an integral of u{T\Kf, decays 
exponentially at positive infinity. Finally, substituting (6.18) into this expression and integrating by parts 
(keeping in mind the above mentioned decay of r(T)), we obtain (1.13). 


32 




REFERENCES 


References 


Acknowledgements 

The work of A. Bogatskiy was supported by the SPbGU grant 11.38.215.2014. T. Claeys was supported by the 
European Research Council under the European Union’s Seventh Eramework Programme (PP/2007/2013)/ 
ERC Grant Agreement n. 307074 and by the Belgian Interuniversity Attraction Pole P07/18. The work of A. 
Its was supported in part by NSE grant DMS-1361856 and the SPbGU grant 11.38.215.2014. 


References 

[1] M. J. Ablowitz and H. Segur. Asymptotic solutions of the Korteweg-deVries equation. Studies inAppl. Math., 
57:1 (1976/77), pp. 13^4 (cit. on p. 3) 

[2] J. Baik, R. Buckingham, and J. DiFranco. Asymptotics of Tracy-Widom distributions and the total integral of a 
Painleve II function. Comm. Math. P/iy^., 280:2 (2008), pp. 463-497. doi: 10.1®®7/sQQ220-0Q8-®433-5 (cit. 
on p. 4) 

[3] J. Baik, R. Buckingham, J. DiFranco, and A. Its. Total integrals of global solutions to Painleve II. Nonlinearity, 
22:5 (2009), pp. 1021-1061. doi: 1®. 1088/0951-7715/22/5/006 (cit. onp. 4) 

[4] M. Bertola. On the location of poles for the Ablowitz-Segur family of solutions to the second Painleve equation. 
Nonlinearity, 25 \A {10\l),pp. 1179-1185. doi: 1®. 1088/0951-7715/25/4/1179 (cit. onp. 3) 

[5] O. Bohigas, j. X. DE Carvalho, and M. P. Pato. Deformations of the Tracy-Widom distribution. Phys. Rev. E 

(3), 79:3 (2009), pp. 031117, 6. doi: 1®. 1103/PhysRevE. 79.031117 (cit. onpp. 9, 10) 

[6] O. Bohigas and M. Pato. Missing levels in correlated spectra. Physics Letters B, 595:1-4 (2004), pp. 171-176. 
doi: http: //dx. doi .org/10.1016/j .physletb. 2004.05.065 (cit. on pp. 9, 10) 

[7] T. Bothner. Transition asymptotics for the painleve ii transcendent (2015). arXiv: 1502.03402 [math-ph] (cit. 
on pp. 4, 10) 

[8] T. Bothner, P. Deift, A. Its, and I. Krasovsky. On the asymptotic behavior of a log gas in the bulk scaling 
limit in the presence of a varying external potential, (2014). arXiv: 1407.291® (cit. on p. 10) 

[9] T. Bothner and A. Its. The nonlinear steepest descent approach to the singular asymptotics of the second Painleve 
transcendent. Phys. D, 241:23-24 (2012), pp. 2204-2225. doi: 10.1016/j .physd. 2012.02.014 (cit. on p. 5) 

[10] Y. Chen and G. Pruessner. Orthogonal polynomials with discontinuous weights. J. Phys. A, 38:12 (2005), 
pp. L191-L198. doi: 10.1088/0305-4470/38/12/L01 (cit. on p. 6) 

[11] P. A. Deift. Orthogonal polynomials and random matrices: a Riemann-Hilbert approach. Vol. 3. Courant Lecture 
Notes in Mathematics. New York University, Courant Institute of Mathematical Sciences, New York; American 
Mathematical Society, Providence, RI, 1999, pp. viii-i-273 (cit. on pp. 6, 12, 25) 

[12] P. A. Deift, A. Its, and I. Krasovsky. Asymptotics of the Airy-kernel determinant. Comm. Math. Phys., 278:3 
(2008), pp. 643-678. doi: 1®. 1007/s®022®-®07-0409-x (cit. on p. 4) 

[13] P. A. Deift and X. Zhou. A steepest descent method for oscillatory Riemann-Hilbert problems. Asymptotics 
fortheMKdV equation. Ann. of Math. (2), 137:2 (1993), pp. 295-368. doi: 10.2307/294654®. arXiv:math/ 
9201261 (cit. onpp. 18,23) 

[14] P. A. Deift, T. Kriecherbauer, K. T.-R. McLaughlin, S. Venakides, and X. Zhou. Strong asymptotics of 
orthogonal polynomials with respect to exponential weights. Comm. Pure Appl. Math., 52:12 (1999), pp. 1491- 
1552. doi: 10.10®2/(SICI) 1097-0312(199912)52 : 12< 1491: :AID-CPA2>3.3.CO;2-R (cit. onp. 22) 


33 



REFERENCES 


References 


[15] P. A. Deift, T. Kriecherbauer, K. T.-R. McLaughlin, S. Venakides, and X. Zhou. Uniform asymptotics for 
polynomials orthogonal with respect to varying exponential weights and applications to universality questions in 
random matrix theory. Comm. Pure Appl. Math., 52:11 (1999), pp. 1335-1425. doi: 10.10Q2/(SICI)1097- 
0312(199911)52:11<1335: : AID-CPA1>3.0 .CO; 2-1 (cit. onp. 6 ) 

[16] P. Deift and D. Gioev. Universality at the edge of the spectmm for unitary, orthogonal, and symplectic ensembles 
of random matrices. Comm. Pure Appl. Math., 60:6 (2007), pp. 867-910. doi: 10.1002/cpa. 20164 (cit. on 

p. 11 ) 

[17] P. A. Deift, A. Its, and I. Krasovsky. Asymptotics of Toeplitz, Hankel, and Toeplitz+Hankel determinants with 
Fisher-Hartwig singularities. Ann. of Math. (2), 174:2 (2011), pp. 1243-1299. doi: 10.4007/annals . 2011. 
174.2.12. arXiv:0905 .0443 [math. FA] (cit. onp. 7) 

[18] F. J. Dyson. “The Coulomb fluid and the flfth Painleve transcendent”. In: Chen Ning Yang. Int. Press, Cambridge, 
MA, 1995, pp. 131-146 (cit. onp. 10) 

[19] A. S. Fokas, a. R. Its, and A. V. Kitaev. The isomonodromy approach to matrix models in 2D quantum gravity. 
Comm. Math. Phys., 147:2 (1992), pp. 395^30 (cit. on pp. 6 , 12) 

[20] A. S. Fokas, A. R. Its, A. A. Kapaev, and V. Y. Novokshenov. Painleve transcendents. Vol. 128. Mathematical 
Surveys and Monographs. The Riemann-Hilbert approach. American Mathematical Society, Providence, RI, 2006, 
pp. xii+553. Doi: 10.1090/surv/128 (cit. onp. 7) 

[21] A. Its and I. Krasovsky. “Hankel determinant and orthogonal polynomials for the Gaussian weight with a jump”. 
In: Integrable systems and random matrices. Vol. 458. Contemp. Math. Amer. Math. Soc., Providence, RI, 2008, 
pp. 215-247. Doi: 10.1090/conm/458/08938. arXiv:0706 .3192 [math. FA] (cit on pp. 2, 4-6, 11, 12, 14,27) 

[22] A. R. Its, A. B. J. Kuijlaars, and J. Ostensson. Asymptotics for a special solution of the thirty fourth Painleve 
equation. Nonlinearity, 22:7 (2009), pp. 1523-1558. doi: 10.1088/0951-7715/22/7/002. arXiv:0811 .3847 
[math.CA] (cit. onpp. 7, 18, 20) 

[23] A. R. Its. “Discrete Painleve Equations and Orthogonal Polynomials”. In: Symmetries and Integrability of 
Difference Equations. Ed. by D. Levi, P. Olver, Z. Thomova, and P. Winternitz. London Mathematical 
Society Lecture Note Series 381. Cambridge University Press, July 2011. Chap. 5, pp. 139-159. doi: 10.1017/ 
CBO9780511997136.007 (cit onp. 5) 

[24] A. R. Its. “Large N asymptotics in random matrices: the Riemann-Hilbert approach”. In: Random matrices, 
random processes and integrable systems. CRM Ser. Math. Phys. Springer, New York, 2011, pp. 351-413. doi: 
10.1007/978-l-4419-9514-8_5 (cit onp. 12) 

[25] K. Johansson. On fluctuations of eigenvalues of random Hermitian matrices. Duke Math. J., 91:1 (1998), 
pp. 151-204. Doi: 10.1215/S0012-7094-98-09108-6 (cit onpp. 2, 6 ) 

[26] A. Kapaev. Global asymptotics of the second Painleve transcendent. Phys. Lett. A, 167:4 (1992), pp. 356-362. 
doi: 10.1016/0375-9601(92)90271-M (cit onpp. 4, 5, 7) 

[27] 1. V. Krasovsky. Correlations of the characteristic polynomials in the Gaussian unitary ensemble or a singular 
Hankel determinant Duke Math. J., 139:3 (2007), pp. 581-619. doi: 10.1215/S0012-7094-07-13936-X. 
arXiv:math-ph/0411016 (cit. onp. 29) 

[28] M. L. Mehta. Random matrices. Third. Vol. 142. Pure and Applied Mathematics (Amsterdam). Else¬ 
vier/Academic Press, Amsterdam, 2004, pp. xviii- 1-688 (cit. on pp. 9, 11) 

[29] M. Plancherel. Sur les valeurs asymptotiques des polynomes d’Hermite Hn{x) - {-Ife^ ^ j. Comment. 
Math. Helv., 1:1 (1929), pp. 227-254. doi: 10.1007/BF01208365 (cit onpp. 5, 6) 

[30] D. Romik. The Surprising Mathematics of Longest Increasing Subsequences. Institute of Mathematical Statistics 
Textbooks. Cambridge University Press, 2015 (cit. on p. 10) 


34 



REFERENCES 


References 


[31] C. A. Tracy and H. Widom. Level-spacing distributions and the Airy kernel. Comm. Math. Phys., 159:1 (1994), 
pp. 151-174 (cit. onp. 3) 

[32] S.-X. Xu, D. Dai, and Y.-Q. Zhao. Critical edge behavior and the Bessel to Airy transition in the singularly 
perturbed Laguerre unitary ensemble. Comm. Math. Phys., 332:3 (2014), pp. 1257-1296. doi: 10. lOO7/sO022O- 
014-2131-9 (cit. onp. 20) 

[33] S.-X. Xu and Y.-Q. Zhao. Painleve XXXIV asymptotics of orthogonal polynomials for the Gaussian weight with 
a jump at the edge. Stud. Appl. Math., 127:1 (2011), pp. 67-105. doi: 10.1111/j . 1467-9590.2010.00512.x 
(cit. onpp. 1-3, 5-7, 13, 14, 16, 18, 24) 


35 



